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We study the effects of low-energy nodal quasiparticles on the classical phase fluctuations in a 
two-dimensional d-wave superconductor. The singularities of the phase-only action at T — > are 
removed in the presence of disorder, which justifies using an extended classical XY-model to describe 
phase fluctuations at low temperatures. 

I. INTRODUCTION 

The spectacular successes of the Bardeen-Cooper-Schricffer (BCS) mean-field theory of superconductivity are based 
on the fact that the so-called Ginzburg-Levanyuk number Giijy), which controls the size of fluctuation effects in D- 
dimensional samples, is very small in bulk conventional materials^ On the other hand, the order parameter fluctuations 
become more pronounced, even dominant, in low-dimensional systems with a small Fermi energy, e.g. in quasi-two- 
dimensional cuprate superconductors^ In particular, the fluctuations of the order parameter phase in the underdoped 
cuprates are enhanced due to a low value of the superfluid density, leading to the large deviations from the BCS 
picture, including the pseudogap phenomenon^ According to Emery and Kivelson^ the Cooper pairs survive in 
underdoped cuprates far above the critical temperature T c , but without global phase coherence, which is destroyed 
by thermal phase fluctuations through the Berezinskii-Kosterlitz-Thouless mechanism^ This idea has been further 
elaborated by many authors, for a review see, e.g., Ref£. As temperature is lowered, a long-range phase coherence sets 
in, but the phase fluctuations continue to play important role in the superconducting state, remaining predominantly 
classical well below T C A& One might expect however that quantum phase fluctuations eventually take over at the 
lowest temperatures. The crossover temperature T c i between the classical and quantum regimes is quite high (of 
the order of T c ) in clean charged systems, but can be significantly reduced in the presence of dissipation. While the 
estimates of T c ; in cuprates obtained by different groups vary considerably, see e.g. Refs£2ii£, here we adopt the view 
that the dissipation is strong enough to make the quantum effects negligible. 

In this brief review we develop an effective long-wavelength theory of the classical phase fluctuations in <i-wave 
superconductors, both with and without elastic disorder. At the Gaussian level, the lowest-order term in the gradient 
energy expansion is simply p s v 2 s /2, where p s is the superfluid mass density and v s is the superfluid velocity. We 
show that the higher-order gradient terms, which contain Vd s , are singular in a clean system at low temperatures 
due to the presence of gap zerosjii and also discuss the effects of disorder on those singularities. As a by-product of 
our theory, we address the question whether using the classical XK-model in <i-wave superconductors can be justified 
from microscopic theory. 

The article is organized as follows. In Sec. |H] the general field-theoretical description of the bosonic excitations 
in superconductors and a phase-only effective action are derived in the clean case. In Sec. IIIII we focus on the case 
of a two-dimensional neutral d-wave superconductor, which is treated in the nodal approximation. The microscopic 
expressions for the energy of fluctuations are compared to the predictions of the classical XY-mode\. In Sec. IIVI we 
derive the d-wave phase-only action in the presence of elastic impurity scattering. Sec. concludes with a discussion 
of our results. 

II. EFFECTIVE FIELD THEORY: CLEAN CASE 

The starting point of our analysis is the tight-binding Hamiltonian 

H = ^Crr'4cr c rV + ^ U r cl a c ra - g ^ Bl r ,B rr , + - ^2{n r - n )V rr >(n r > - n ), (1) 

rr' r (rr') rr' 

where r label the sites of a tetragonal lattice. The hopping amplitude t rr i and the chemical potential fj, are combined 
into the band-dispersion matrix t; rr i = —t rr i — fiS rr ' , which is real and symmetric in the absence of external magnetic 
field. The second term describes impurity scattering. The third term is the BCS interaction in the <i-wave channel, 
g > is the coupling constant, and the operator B rr ' — {c r 'ic r } — c r '-\c r i)/V2 destroys a singlet pair of electrons at 
the nearest-neighbor sites (rr') in the xy plane. The last term describes the repulsive interaction between electrons, 
n r = c rr7 c r a is the particle number density, no ~ (n r ) is the average number of particles per site (which is equal to 
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the ionic background density, thus ensuring the overall charge neutrality of the system), and V rr > is the interaction 
matrix. 

Let us first look into the clean case. Setting U r = in Eq. l|T|l. the partition function can be written as a functional 
integral over Grassmann fields c rcr {r) and c ro -(r): 

Z = Tre- (3H = J VcVce- s[ec] , (2) 

where S — dr ^2 r c ra d T c ra + H(t)], f3 = l/T (in our units ks = h = 1). Using the Hubbard-Stratonovich 
transformation to decouple the interaction terms, we end up with the representation of Z as a functional integral over 
c, c and two bosonic fields: a complex field A rr i (r) , which describes the superconducting order parameter fluctuations 
and is non-zero only on the bonds between the nearest neighbors, and a real scalar potential field ip r (r). The fermionic 
part of the action then becomes S = Tr(C't/ _1 C), where 

_ ( c rT 



C r = I J , C r = ( C rT C rl ) 

are the Nambu spinors, and Q^ 1 is the inverse Green's operator: 

fi-l/V „' ~>\ — xfr _ ( d rr i[— d T — i<p r (T)] — £ rr > —A rr i(T) \ , , 

y (r,T,r,r)-d[r r)^ -A* rr ,(r) 8 rr , \-d T + i Vr {r)\ + £„, ) ' [6 > 

We use the notation "Tr" for the full operator trace with respect to both the space-time coordinates and the Nambu 
matrix indices, reserving "tr" for a 2 x 2 matrix trace in the Nambu space. Integrating out the fermionic fields, we 
obtain 

Z = [vA*VAVlp e -Se//[A*,A^ 



with the effective action 

S e ff = -Trliif? -1 + J dr l-^\A rr ,\ 2 + ^^ip r V~Jtp r > - in ^2(p r \ . (5) 

q \ rr' rr' r J 

The mean-field BCS theory corresponds to a stationary and uniform saddle point of the effective action lJ5]l. which 
is found from the equations 6S e ff/SA* = 5S e ff/5tp = 0. The solution describing d-wave pairing is given by ipo,r = 
and 

Ao^ = | + ^°' * r 'r r t a *> (6) 
1 — A , if r = r ± ay. K ' 

In the momentum representation, = Ao(T)4>k, where <\>k = 2(cosfc £! ;a — cosk y a) is the <i-wave symmetry factor. 
The temperature dependence of the gap amplitude Ao is determined by the standard BCS self-consistency equation, 
generalized to the case of an anisotropic order parameter^ 

Inverting the operator ©, we find the mean-field matrix Green's function: 

r /, s_( G (k,uj n ) -F a (k,u n ) \ _ iwuTp + $ fc T 3 + A fc n 

' UJn) ~ \-F (k,u; n ) -G (-k,-u> n ) ) - a;2 + £2 +A 2 ■ ^ 

Here u n — (2n+l)nT is the fermionic Matsubara frequency, Go and Fq are the usual normal and anomalous Gor'kov's 
functions of the superconductor^ £ fc = £_ fc is the band dispersion of free electrons, and are Pauli matrices in the 
Nambu space. The Green's function J7J) determines the single-particle properties in the mean-field approximation. 
In particular, after its analytical continuation to the real frequency axis, iu> n — > to + iO, one obtains the energies of 
elementary fermionic excitations, or the Bogoliubov quasiparticles: 



E k = ^l + Al (8) 

For the ei-wave order parameter, the gap in the excitation energy vanishes along the diagonals of the Brillouin zone. 
These zeros, or the gap "nodes" , are responsible for many peculiar thermodynamical and transport properties in the 
superconducting state on the mean- field level. — Below we show that the gap nodes also have dramatic effects on the 
long-wavelength behavior of phase fluctuations. 
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A. Fluctuations 



Deviations from the mean-field solution can be represented in terms of the amplitude and phase fluctuations: 
A rP ,(r) = [Ao, rr / + ( 5A rW (r)]e te -'( r ), where 6 A and are real. We neglect the amplitude fluctuations because they 
are gapped, see e.g. Ref^, and therefore make a negligible contribution at low temperatures. Since the number of 
bonds in a square lattice is twice the number of sites, one needs two on-site phase fields r (r) and d r {r) to describe 
the phase degrees of freedom. One possible parametrization is 



©r 



r , if r' = r + ax 
+ 6 r , if r' = r + ay. 



(9) 



The fluctuations of 0, which describe a change in the symmetry of the order parameter from a pure d-wave to a 
d + is-wave, can be neglected^ 

If r changes slowly over the lattice constant, then one can make the replacement 



A , 



An ,P J ( e -+^')/2 



(10) 



(recall that r and r' are nearest neighbors). The next step is to perform a gauge transformation to make the 
off-diagonal elements of Q~ x in Eq. real: 



U\r,T)Q- x {r,T-y ,T')U{r' ,t') = g-\r,T;r' ,r'), 



(11) 



where U(r,r) — exp[zT30 r (r)/2]. This transformation leaves the operator trace in the effective action (JjJJ invariant. 
Although the order parameter I|1U|) is no longer invariant under local phase rotations r — > r + 2tt, our results are not 
affected since we consider only small fluctuations of the phase in the low-temperature limit, where the contribution 
of vortices can be safely neglected. 

The operator Q" 1 can be represented in the form Q^ 1 = Q^ 1 — S, where Go ls the mean-field matrix Green's 
function, whose Fourier transform is given by Eq. |7|. and 

ld0 r (T) 



E(r, r;r',r') = S(t - t')t 3 \ id. 



2 dr 



-iT 3 [e T .(T)-e 7 .,(r)]/2 _ 1 



(12) 



is the self-energy correction due to fluctuations. At slow temporal and spatial variations of and small tp, one can 
expand the effective action |JSJ in powers of S, keeping only the two lowest orders in the expansion, with the following 
result: 



/9 
S eff [9,cp] =S + Tr(g S) + iT>(a £goE) + 0(S 3 ) + i J dr ^ ^(r)V;> r ,(r) - ^n J dr^^(r), 



(13) 



where 



So = -Tr In G^ 1 + 0~ ^ | A . rr ' | 2 = PS 



(14) 



is the saddle-point action, So is the total mean-field energy of the superconductor, and Go is the saddle-point Green's 
function, see Eq. J7J. For non-interacting fluctuations we keep only the terms of the first and second order in S. 
Calculating the traces in Eq. (|13|) . we obtain the Gaussian action 



S e ff[0,(p] = S Q + Si + S 2 , 



(15) 



where 



dO r (T) 

dr 



ITTHo 



E 



r {(3) - 0r(0) 

2tt 



(16) 



is the topological term containing the phase winding numbers, and 



Z Q 



L W (Q)\^(Q)\ 2 + L ev (Q)0*(QMQ) + -L gB (Q)\0(Q)\' 



(17) 
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Here we use the shorthand notations Q = (q, v m ) and 



£(■••) = t££(...), 



where v m = 2mirT is the bosonic Matsubara frequency and the momentum summation goes over the first Brillouin 
zone. The lattice Fourier transforms of the fields are defined by the usual expressions: 9 r (r) — Af^ 1 ^ 2 Yl q ^0?' T ) elqT 
etc, where M is the number of lattice sites. Since both ip and 9 are real, they satisfy 9*(Q) — 9(—Q), <p*(Q) — <p(—Q). 
The coefficients in Eq. (|17|l are given by 



L vv (Q) = V- 1 (q)+Il (Q), 
LeipiQ) = iv m n (Q) + q t n\(Q), 
Lee{Q) = t^n (Q) - 2iv m q$? 1 {Q) + q iqj Il%(Q). 
Here V (q) is the Fourier transform of the interaction matrix V rr i , and 

n (Q) = - £>[S (# + Q)t 3 Qo(K)t 3 ], 



K 



KiQ) = J2 v MGo(K + Q)ng (K)T ], 

K 

n? (Q) = £ ™^IGo(K)t 3 ] + m tx[G (K + Q)t G (K)t ] 

K 

In these expressions, K — (fc,cj n ), 



(18) 

(19) 
(20) 
(21) 



A" 



£(•••) = ?T^ £(...) ^T^aJ 

K n k n 



(2ir) D 



(•••), 



m 7j 1 (^) = d 2 £,k/dkidkj is the inverse effective mass tensor, v(k) — d^k/dk is the quasiparticle band velocity, and £2 
is the unit cell volume (to simplify the notations, below we set £2 = 1). Eqs. (|18|) are obtained in the limit of small 
q from more general expressions, using the gradient expansion £fc+q = £fc + vq + (l/2)m~j 1 qiqj + 0(q 3 ). Integrating 
out the field ip, we finally arrive at the phase-only effective action: 



Seff[9} = So + S 1 + lYl 



Q 



L ee {Q) 



L 2 8v (Q) 



L VV (Q) 



(22) 



In this article, we focus on the case of classical phase fluctuations and neglect all interactions other than those 
responsible for the Cooper pairing, which corresponds to neglecting the r-dependence of 9 and setting V(q) = 0. 
Then the topological term vanishes and only the iSf contribution survives in the third term in Eq. <|22[) . so that the 
effective action becomes 



S ef f[0\=l3£o + I3£[0], 



(23) 



where £ is the energy of fluctuations in the Gaussian approximation. Calculating the Matsubara sums in Eq. 121|) 
and introducing the superfluid velocity v s = (l/2m)V#, where m is the electron mass, we obtain 



£ = 5 £%(g)<;(gK,j(£?), 



(24) 



with the kernel 



where 



K. 



(0) 



2m 2 T]T J 



d D k 



(27T) 



Jdj(q) = m 2 H^( q ,0) = JC[f +lij(q), 



m i 1 {k)Go(k,uj n ), 



D i] 



d D k 



Vi(k)vj{k) 



■ g) tanh %*+ tanh # + ^ tanh - tanh # 



E, 



k+q 



+ E k 



Ek+q — E, 



k 



(25) 

(26) 
(27) 
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and 

C ± (k, q ) = Ul± ^ + F AkAk+ A (28) 

are the coherence factors. 

An important characteristic of the superconductor is the superfiuid density tensor, which is defined as 

p a #(T) = Kij(0)- (29) 

Its temperature dependence can be easily found in two limiting cases. In the normal state, A& = 0, and one can use 
the identity dGo/dk — vGq in Eq. (|21(l to obtain p s ,ij{T > T c ) — 0. On the other hand, at zero temperature we have 

p s ,y(0) = K,f), since 2^(0) = at T = 0. 

The expressions (I25I26I27[) are valid for arbitrary band structure and gap symmetry In a Galilean-invariant 
system, i.e. for = k 2 /2m — p, the tensor l|2t)[l takes a particularly simple form: fcyy = po%, where 
po = 2mT ^2 n J k Go(k, uj n ) is the mass density of electrons. Therefore, the superfiuid density tensor at zero tem- 
perature is p s ,ij(0) — Po$ij, i-e. all electrons are superconducting. In general, there is no such simple relation in a 
crystal. 



III. TWO-DIMENSIONAL CASE 



In this section, we apply the general theory developed above to a two-dimensional <i-wave superconductor. In this 
case, the low-energy physics at T — ► can be conveniently described using the so-called "nodal approximation" 
which takes advantage of the fact that the excitation energy © for the <i-wave order parameter can be linearized in 
the vicinity of the four gap nodes located at k n = k F k n (n = 1, 2, 3, 4) on the Fermi surface. Here 

£ x + y - -x + y -x-y } x~y 
fe i = — 7^~' fe 2 = j= — , k 3 = ■= — , fe 4 = — -=-. 

V2 y/2 V2 V2 

For the nodal quasiparticles in the vicinity of the nth node we have k = k n + Sk, and 



£fc = v F 5k ± , A fe = v A 6k\\, E k = Jv 2 F Sk 2 j_ +vlSkf,, (30) 



II 

where Sk± and 8ku are the momentum components perpendicular and parallel to the Fermi surface, v F is the Fermi 
velocity at the nodes, and va = \dAk/dk\\\ is the slope of the superconducting gap function near the nodes. Thus, the 
excitation spectrum near the gap nodes is described by an anisotropic Dirac cone. The anisotropy ratio u_f/i>a is an 
important characteristic of the high-T c cuprates, which depends on the material and the doping level, e.g. v f /va — 14 
and 19 in the optimally doped YBCO and Bi-2212, respectively^ 

In the nodal approximation, the momentum integration over the whole Brillouin zone is replaced by a sum of four 
integrals over the small regions in fe-space around the nodes: 



(I 



In the last integral, we changed to the polar coordinates: VpSkx — ecosa, «A^fc|| = esina, and E k = e. The 
ultraviolet cutoff e max ~ Ao is introduced to make sure that the area of the integration region is equal to the area of 
the original Brillouin zone^ In most calculations in this article this cutoff can be extended to infinity. 



A. Classical phase fluctuations 

We can now calculate the energy of the classical phase fluctuations, see Eq. I|24() . The nodal approximation cannot 
be applied to the momentum integral in Eq. i|26|) because it contains contributions from all electrons, including those 
far from the Fermi surface. Assuming that the band dispersion can be treated in the effective mass approximation, 
which amounts to the replacement — > (l/m*)<5y , one obtains 

= (32) 
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where po is the average mass density of electrons. 

In contrast, the second term in the kernel IC can be calculated in the nodal approximation. Using Eqs. (|3UI31[) . we 
have 



Tij{q) 



2tt v a ^ 



(33) 



where 



Here 



Si(q) = S 3 (q) = Ts , S 2 (q) = S^q) = Ts (|) 



7i,2(q) = -^Jv 2 F {q x ± gy ) 2 + «|(g x T 9y) 2 
are the energies of the nodal quasiparticles with 5k = q, and the scaling function s(x) is defined by an integral: 

OO 2-7T 



(34) 



s(x) 



dyy 



da 
2^ 



[f+(x,y,a) + f-(x,y,a)), 



(35) 



where 



1 ± 



y + xcosa \ tanh(\/a; 2 + y 2 + 2xy cos a/2) =p tanh(y/2) 
2xy cos a y \/a; 2 + y 2 + 2xy cos a =F y 



\Jx 2 + y 2 



Note that the cutoff energy e max has been replaced by infinity, due to the rapid convergence of the integrals. One can 
show that the function s{x) has the following asymptotics: 



s(x) 



21n2 + ^ 

7TX 



at x 
at x 





00. 



(36) 



After the summation over the four nodes in Eq. 1)33(1 . one finally obtains 



r („\ - m n a ™ 2 Vf t ( 



(37) 



where 



*±(g) 



"71(g)" 


± s 


'72(g)" 


T 




T 



The expression (|37|l is exact in the nodal approximation, i.e. for the conical quasiparticle spectrum. In terms of q, the 
applicability region of the nodal approximation is 71,2(9) Ao. At higher energies of quasiparticles, the deviations 
of the spectrum from the linearized form (|30|l should be taken into account. 

We would like to note that in the nodal approximation, U\(q,Q) = and therefore Lg ip (q,0) — 0, see Eqs. <|18ll . 
(I20|l . This means that the classical fluctuation energy in two dimensions has the form (|24l) with the kernel H37[). even 
if the Coulomb interaction is taken into account. 

At T — 0, using the large- a; asymptotics of s(x) in Eq. I)36|) . the kernel takes the form 



m 2 v F ( 7l (q) + 72 (q) 7l (q) - 72 (q) 



13 16 v A V 7i(g) ~ 72(g) 7i(g) + 72(g) 



(38) 



Setting q = here, we find the superfluid density: p s (Q) = (jn/m*)po. We also see that the kernel is a non-analytical 
function of q, which means that no gradient expansion of the energy exists. 

At finite temperatures and small q, such that 71.2(g) -C T, the small- x asymptotics of s(x) yields 



, ..1 2 In 2 v F o™ , - 

^i(g) = — Po — m T ) 6 i4 - — — - 

m* it va 



m 2 V F 1 / (4 + v ^)( q l + q2) 2(v 2 F - v 2 A )q x q y 



itvaT\ 2(v 2 F - v 2 A )q x q y {v 2 F + v\){q 2 x + q 2 ) 



(39) 
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The first term describes the depletion of the superfluid density due to the thermal excitation of quasiparticles: 

, > , - 2 In 2 vp o , 
p s (T) = p s {0) ^m 2 T, 40 

TT VA 

see also Refs»i£*iii, which explains a linear in T increase of the magnetic penetration depth X(T) at low 
temperaturesj22i2i observed in high-T c cupratesi^ 

The quadratic q-dependence of the expression (|39J) implies that the kernel JCij (R) in real space is proportional to 
exp(— with the length £ given by 

fm= C^LV^^V /2 & , (41) 



\A8TrvAPsT J \T 

where £o = vf/^ttT c is the BCS coherence length (we assumed that vf 3> va)- This behavior is similar to that of the 
electromagnetic response function in conventional s-wave superconductors, see e.g. Ref^. An important difference 
however is that the characteristic length £ is temperature-dependent: £(T) ~ 7 1-1 / 2 at T — > 0. It is because of the 
divergence of £ that the gradient expansion of the classical energy of fluctuations breaks down at T — ► 0. Note also 
that the length £ is different from other characteristic lengths discussed in the literature: £o ~ the coherence length, 
or the correlation length of the gap amplitude fluctuations, which remains constant at T — > 0, and £ pa ir ~~ the size of 
a Cooper pair, which is infinite in the d-wave case^ In a conventional s-wave superconductor, all three lengths are 
of the same order. 

The physical interpretation of our findings is the same as that of the non-local Meissner effect^ since the gap 
function A*, has nodes on the Fermi surface, then the anisotropic coherence length t\p/|Afc| exceeds the London 
penetration depth Xl close to the nodes, and the local electrodynamics breaks down. 

One can expect that the nodal quasiparticles also affect the non-Gaussian terms in the effective action As 
shown in the Appendix, indeed no expansion of the classical fluctuation energy in powers of v s exists at T = 0. 



B. Failure of the classical XF-model 



The effects of the phase fluctuations in superconductors are most often studied using either the classical or the 
quantum versions of the XY-mode\. The energy of the classical AF-model in the absence of external fields has the 
form 

SxY=^2jRR'[l-ooa(eR-e w )], (42) 

RR' 

where Or is the phase of the order parameter at site R of a coarse-grained square lattice, whose lattice spacing d is 
of the order of the superconducting correlation length £. The summation goes over all bonds in the lattice, and the 
coupling constants Jrr> — J(p), where p = R — R' , are called the phase stiffness coefficients. While in the Gaussian 
approximation, see below, the coupling constants are temperature-independent, the interaction of fluctuations leads 
to a thermal renormalization of the Js, and eventually to a phase transition into the disordered state. 

The XK-model is believed to provide a correct description of any system with broken U(l) symmetry if the 
amplitude fluctuations of the order parameter are negligible. Typical examples are classical Heisenberg magnets, 
superfluids and superconductors. The experimental systems to which the lattice model (|42H has been applied include 
granular superconductors and fabricated arrays of Josephson junctions, see, e.g., Refj2^ and the references therein. 
In those cases d is given by the distance between grains. Although the simplicity of the AY-model is physically 
appealing, its rigorous microscopic derivation for homogeneous high-T c superconductors does not exist. The usual way 
of justification, see, e.g., Ref^, involves expanding the cosine in Eq. I|42(l and matching the expansion coefficients with 
those in the Gaussian phase-only action. That the microscopic theory fails to reproduce the quantum generalization 
of the XF-model has already been noticed in Refi^&. Here we show, following Refiii, that even for the classical phase 
fluctuations in a <i-wave superconductor at low temperatures, the long- wavelength limit of the microscopic theory is 
not consistent with Eq. 142|) . 

For slow variations of the phase, the energy (|42|l takes a Gaussian form: 

RR' q p 
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In terms of the superfluid velocity, we have 

£ *y = l'E lc v r bKM v >Aq)- ( 43 ) 
q 

The kernel here has a well-defined Taylor expansion in powers of q: 

£* Y (q) = p*Jj + Kmwi + o(<? 4 ), (44) 

where 

pXJ 3 =2m 2 Y J J{p)p l P 3 (45) 
p 

is the superfluid mass density tensor (for example, if the only non-zero coupling is between the nearest-neighbor sites, 
then pfji = 8m 2 d 2 JSij), and 

2 

Mj,kl = — g- 22J(p)piPjPkPl- (46) 
p 

Comparing Eqs. (|44|l and (|38|) . we see that the effective long- wavelength theory of the classical phase fluctuations 
at T = does not have the form of the XF-model, since the momentum dependence of the two energies is clearly 
different. At T > 0, although the expression (|39[) is quadratic in q, the coefficients diverge as T — > 0, which is not the 
case for A^jy above. Thus, the microscopic theory fails to reproduce the long-wavelength structure of the classical 
AY"-modcl in a clean d-wave superconductor. 



IV. DISORDERED CASE 



In the presence of impurities, a full effective field theory for the disordered interacting system described by the 
Hamiltonian JTJ would include the fluctuations of the order parameter and of the scalar potential coupled with the 
disorder- induced soft modes (the diffusons and the Cooperons). Such theories, usually having the form of a non-linear 
cr-model, have been developed, see, e.g. Ref^ and the references therein, to study the effects that are beyond the 
scope of the present work, for instance the suppression of T c due to the interplay of disorder and interactions in s-wave 
superconductors. Our goal here is to check if the elastic impurity scattering removes the divergencies in the gradient 
expansion of the classical phase-only action discussed above. The disorder is treated essentially in the saddle-point 
approximation and the Coulomb interaction is neglected. 

As a bookkeeping device to obtain an effective action for the order parameter fluctuations, we use the replica trick: 
(InZ) = \im n ^o((Z n ) — l)/n (the angular brackets denote averaging with respect to disorder). From Eq. Q we have 



/n 
Y[Vc a Vc a e- sl5 ' c] , (47) 



where S = 5 + S int , 



/a 

5 ° = E / dr (E 4«rC + E + E ^cc) , 



(48) 



o 

P 



Sint = -9 E / dT E S rr' B rr< ■ (49) 
a q (rr'> 

The impurity potential here is assumed to be Gaussian-distributed, with zero mean and the correlator 

(U r U r i) — — S rr i , 
zttNft 
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where Np is the density of states at the Fermi level, and r is the electron mean-free time due to elastic scattering. 
The next step is to use an incomplete Hubbard-Stratonovich transformation to decouple only the interaction terms 
Si n t in each replica, before disorder averaging^ Proceeding as in Sec. [H] we have 

/n 
Y[VA a '*VA a e- s '"^ A "' A \ (50) 



with the effective action 



p 

S eff = -]n(Detg- 1 ) + -T t f dr^Kr'l 2 , 



(51) 



instead of Eq. J5J. Here 



and "Det" stands for the full operator determinant with respect to the space-time coordinates and the Nambu and 
the replica indices. The disorder averaging in the first term in S e ff generates effective coupling between different 
replicas, resulting in a non-linear bosonic field theory. 

While the saddle point of the effective action l|51|l has the same structure as in the clean case: Ag , (r) = A rr i , 
see Eq. the temperature dependence of the gap amplitude is different, in particular, the critical temperature is 
suppressed by impurities^ The mean-field Green's function is the unity matrix in the replica space: = Sa^Q^ 1 , 

where 

n-lf„ \ — ( &rr'{iu n — U r ) — £,rr' _ ^0,rr' \ /ro\ 

y [r, r , cj n) - y ^ ^ ^ ( + u j + ^ J W 

in the Matsubara frequency representation. The disorder averaging can be done using the standard diagram 
technique^ In the so-called self-consistent Born approximation, which assumes a sufficiently weak impurity scat- 
tering and also neglects the diagrams with crossed impurity lines, the average mean-field Green's function has the 
form 

<«>(*, «»)> = a2 + g + A | . (54) 



where w n satisfies the equation iu n — iuo n + (1 / 2t) ((iuj n / \/oj 2 + A|))fc>i- 

Assuming replica-symmetric phase fluctuations, the order parameter can be written in the form A° r , (r) = 
Ao j rr'e 4 ' 6 ' T ' (T ' )+6 ' 7, ' ( - T ' ) ]/ 2 , see Eq. Performing the gauge transformation (| 1 1 11 in each replica, we have Q~ b = 

dabiGo 1 — S), where the self-energy operator has the form Ijl2|l with f r (j) — 0. Therefore, 

^(Detg- 1 ) = ln(Det CT 1 ) = ^(Trln^ 1 ) + n(Trln(l - g E)) +0(n 2 ), 

in the limit n — > 0. Substituting this expansion in Eq. I|50|) . we see that the replica index can be omitted, and the 
effective action (TH51) gets replaced by its disorder average: 

S eff [0] = (3£ - (Trln(l - O S)), (55) 

where £q is the mean-field energy of a disordered superconductor. Considering only the static fluctuations and 
expanding the operator trace in powers of the phase gradients, we obtain the effective action in the form (|23[l . where 
the energy of fluctuations is now given by 



£ = 5£<M9)Ki(9Kj(9), (56) 

with 



2 
q 



(lC ij (q)) = (]C^) + {l ij (q))- (57) 
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Before proceeding with the calculation of the disorder averages, we would like to note that Eq. I|56|) could also be 
derived using a less formal approach, without introducing replicas. Assuming that electrons move in the presence of 
the random potential and a given order parameter field A rr / = /S.Q irr 'e l ^ e7,+Br '^ 2 , one can define the energy functional 
of phase fluctuations: 

£[e]=-~{hxZ\e]) + ±{\nZ[Q]), (58) 

where Z[9] = Det (Gq 1 — S) is the partition function. Expanding lnZ[#] in powers of the phase gradients, followed 
by averaging each expansion coefficient with respect to disorder, we arrive at Eq. 15611 . 
The disorder averaging of the first term in Eq. 1)57(1 is straightforward: 



(4 0) ) = 2m 2 Tj2 J ^pm-- 1 (fc)(G (fc,^)> = ^rPo%, (59) 



where we used the effective mass approximation, as in Eq. (|32[) . 

The average of the product of two Green's functions in the second term includes the impurity vertex corrections. 
However, since we are interested in the behavior of the kernel in the long-wavelength limit q — > in a sufficiently 
clean system (the precise criterion will be discussed below), the vertex corrections are negligible. Thus 



/d D k 
^^ D Vi(k)vj(k) tr[{g (k + q,uj n ))T (g (k,uj n ))To} 

In order to calculate the Matsubara sum we use the spectral representation 



(60) 



, . de Im(X?(fc, e) . „ 

g (k,u n )= / y ° K ' ' , (6i) 



where the retarded matrix Green's function is obtained from Eq. I)54|) by analytical continuation: 



and t(e) — {iLu n )\ioj n ^e+io satisfies the equation t = e + (i / '2t) ((t / \J t 2 — A|))fe. Instead of finding the exact energy 
dependence of t(e) for the d-wave gap, we use below a simplified expression 

t( e )=e + iT, (63) 

which captures the essential qualitative effects of the impurity broadening of the single-electron states, with T being 
an energy-independent effective scattering rate. The exact solution shows that the scattering rate indeed approaches 
a constant value at e — ► 0, i.e. in the so-called universal limits 
Inserting the representation (|61Jl in Eq. (|60JI . one obtains 

oo oc 

. , .. m 2 f , r tanh^-tanh^ f d D k . , 
<Z y (q))=-— J de, J de 2 J __ Vi{k)vj{k) 

— oo — OO 

x {C_(fc, q)[d+{k + q, ei)d_(fe, e 2 ) + d-(k + q, e x )d + (k, e 2 )] 

+C+(fe, q)[d+(k + q, e 1 )d + (k, e 2 ) + d_{k + q, ex)d_(fe, e 2 )]}, (64) 

where 

i r 



d±(k, e) = 



7T(e±E k ) 2 + T 2 - 



and C± are the coherence factors (|28|) . The expressions l|27ll are recovered from Eq. (|64|l in the clean limit T — > 0, 
when d±(k, e) — » S(e ± Ek). 

We focus on the case of zero temperature, when the kernel is a non-analytical function of q in the absence of 
impurities, see Eq. (|38|) . At T = 0, the integrals over ei j2 can be calculated, giving 



2m 2 f d D k 



(lij(q)) = — / j^-jD v *(k)vj(k) 



arctan ^±3. _|_ arc tan ^ arctan ^±3. _ arc tan ^ 



C_(fc, g) =± — ^ + C+(fc, g) 



-E'fe+g + E k E k+q — E k 



(65) 
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Comparing this to Eq. I|27|l. we see that the energy of classical phase fluctuations in the disordered case at 
zero temperature has exactly the same form as in the clean case at a finite temperature, if one formally replaces 
tanh(-E/2T) — > (2/w) arctan(i?/T). Therefore one can expect that the disorder will affect the phase fluctuations in 
the same way as temperature does, i.e. the singularities of the effective action will be washed out. 

To check this conclusion quantitatively, let us evaluate the long- wavelength asymptotics of (Tij(q)) in the nodal 
approximation. Using Eqs. (|30I31(I . we have 

m 2 4 

(2ij(q)) = V k n ,ikn,jS n (q), (66) 

Z-K V A , 

where 

Sx(q) = S 3 (q) = TS , S 2 (q) = S±(q) = TS (^) , 

and 71,2(9) are given by Eq. (|34f> . The function s(x) is defined by 

s(x) — I ~z / dyy[f+{x,y,a) + f-(x,y,a)], (67) 



2tt 



where y max = e max /T, and 



~ if ty + xcosa \ arctan(-\/a; 2 + y 2 + 2xy cos a) =F arctan(y) 

J± = - 1 ± 



7T 



\/ x 2 + y 2 + 2xy cos a J \/ x 2 + y 2 + 2xy cos a =F y 



In contrast to the clean case, see Eq. ().'{ot , the energy cutoff here cannot be extended to infinity, due to the logarithmic 
divergence of the integrals. 

The clean-case expression is recovered from Eq. (|fiTjl at x 3> 1, i.e. at T <C 71,2(9) <C fmai- In the opposite 
limit x -C 1, one finds 

2 1 

s(x) = - lny maa; + — a; 2 , (68) 

7T D7T 

which is valid at y ma x 3* 1- Inserting this into Eq. I|tjtj|) . w< 3 finally obtain 

/r ^ - ( m n 2 VF m 2 r^ e,nax \t m2 v f ^ ( {v 2 F +v 2 A ){q 2 x +q 2 y ) 2{v 2 F -v A )q x q y \ 
\m* 7H i; A r / 12tt z v a r \ ^{v F ~ v A )q x q y (v F + v A )(q x + q y ) J „ 

The first term describes the depletion of the superfluid density at zero temperature due to the impurity scattering: 

2 v F 2 A 
— — m Tin — 

7T Z 1 



Ps {T = 0, T) = p s (T = 0, 0) - -2 — m 2 rin (70) 



see also Ref^i (here we used e max ~ Ao). Upon increasing T, the superfluid density decreases and eventually vanishes 
at some critical disorder strength. One cannot reach the critical point using the expression (|70l) . since its validity is 
limited to the case of weak impurity scattering T <C Aq. 

The second term in Eq. l|69fl shows that, in contrast to the clean case, the energy is a non-singular function of 
momentum, even at T = 0. Comparing it to the expression (|44|l . we see that the microscopic theory, at least at 
the Gaussian level, has the same long-wavelength behavior as the classical Xy-model. The non-zero elements of the 
tensor A are given by 



_ _ m 2 1 , 2 2 \ 

^-xx,xx — ^-xx,yy — ~TZ o T^(^_F ^A/ 
lZ7T VA i 

_ m 2 v F 1 , 2 

^xy.xy — 10 _2.. t^\ V F~ v A)i 



"' "' y 12tt 2 v a T' 



with other elements obtained by symmetry. One can see that keeping only the nearest-neighbor phase stiffness 
coefficient in Eq. (|46|l is not sufficient to reproduce the tensor structure (|T1|) . therefore one has to consider an 
extended XY"-model, in which the couplings between next-nearest neighbors etc, are also taken into account. This 
was first noticed in Ref^S for an s-wave superconductor. 
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V. CONCLUSIONS 



To summarize, we have shown that the effective field theory for classical phase fluctuations in a clean rf-wave 
superconductor suffers from singularities which make the gradient expansion of the fluctuation energy impossible. 
This means, in particular, that the physics of classical phase fluctuations cannot be described by the XY-mode\ at 
low temperatures. In the presence of disorder, a well-defined gradient expansion of the Gaussian phase-only action is 
restored, and has the same form as that of an extended classical XF-model. 
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APPENDIX A: CONDENSATE ENERGY AND NON-LINEAR MEISSNER EFFECT 



In this Appendix we calculate the higher-order terms in the expansion of the effective action 

S eJ j = S - TrlnCT 1 + Trln£ ~\ (Al) 

in powers of the phase gradients. We consider only the limit of uniformly-moving condensate, when the superfluid 
velocity v s is constant, so that 9 r — 6 r > = 2mv s (r — r'). In the absence of the scalar potential, the self-energy l|12|l 
becomes translationally invariant: 

E(r, r; r', r') = 6(r - r'^w [e-^m^r-r') _ x j _ 

The gauge-transformed Green's function 1)11(1 becomes diagonal in the momentum space: 

0-\k,o>n) ~t k t-), (A2) 



y -Afe itv n + Cfc / 
where £. = ^± m „ ( , which allows one to calculate the operator traces: 

Trln^- 1 =lnDctg- 1 =/3VT^ J ^5 lndet^-^fe, w„) 

where V is the system volume and "det" denotes a 2 x 2 matrix determinant in the electron-hole space. Inserting this 
in Eq. (HJ, we obtain S e ff — Sq + (3£, where 

£ = -VTY [ In - + %j (A3) 

has the meaning of the kinetic energy of uniformly moving condensate, and 



are the quasiparticle energies affected by the superflow. Using the identity 



E, ito n — E „ , cosh 
In — = Tin 



2T 



iuj n - E cosh ' 



the energy density can be written in the form 



£ f d D k , cosh %t cos h E* 



r- f l^ h ~^» 2T <A5) 
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At T = 0, this becomes 

£ 1 f d D k 



V 2 / (2tt 



(\E h ,+ \ + \E k ,+ \-2E h ). (A6) 



We focus now on the case of a two-dimensional d-wave superconductor. Assuming for simplicity a Galilean-invariant 
case, characterized by the parabolic dispersion = k 2 /2m — /j,, with the effective mass equal to the bare electron 
mass to, we have 



E k ,± = yj (a + C) 2 + A| ± kv s = £ c , fc ± fcu s , 

where fcu s is the so-called Doppler shift of the quasiparticle energy in the presence of moving condensate, and 
£ = mv 2 /2. The expression (|A6() can then be written as £ /V = {£/V) reg + {£ /V)dopp- The first term, 

^ - fd2k {E c , k -E h ) = - Cf^^ + 0(e) = ^ + { v% (A7) 



Vj reg J (2^) 2 v ^ *' *J (2tt)2 E k 

has a well-defined Taylor expansion in powers of v s . Note that the nodal approximation cannot be used here because 
of the contributions from the regions far from the Fermi surface, see Sec. 1111 Al In contrast, the second term, 

\E c , k +r ) \ + \E C!k -r ] \-2E c , k \, (A8) 



( L ) = 




f d 2 k 


\V/ Dopp 




(27T) 2 



can be calculated in the nodal approximation. Using Eqs. (|30I31|I . we obtain 

V J Dopp 12nv F v A ^ s| 
which is a non-analytical function of v s . Putting all pieces together, we arrive at the following result: 

£ = po^ _ rn?^_ + , + _ + 0{<y {M) 

V l 12\/2irvA 

We see that the nodal quasiparticles make the kinetic energy of the condensate a non-analytical function of the 
superfluid velocity^ This singular behavior is closely related to the so-called non- linear Meissner effect^ the screening 
supercurrent acts as a pair-breaker in d-wave superconductors, creating a finite density of normal excitations even 
at T = 0. The depletion of the supercurrent by those excitations leads to a non-analytical dependence of the 
electromagnetic response functions on v s and therefore the external magnetic field. 
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